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We apply results derived in other contexts for the spectrum of the Laplace 

operator in curved geometries to the study of an ideal polymer chain confined 

to a spherical annulus in arbitrary space dimension D and conclude that the 

free energy compared to its value for an uncurved box of the same thickness 
• 5-h , 

and volume, is lower when D < 3, stays the same when D = 3, and is higher 
when D > 3. Thus confining an ideal polymer chain to a cylindrical shell, 
lowers the effective bending elasticity of the walls, and might induce spon- 
taneous symmetry breaking, i.e. bending. (Actually, the above mentioned 
results show that any shell in D = 3 induces this effect, except for a spherical 
shell). We compute the contribution of this effect to the bending rigidities in 
the Helfrich free energy expression. 
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I. INTRODUCTION 



The quantum mechanics of particles confined to surfaces or surface layers embedded 
in a higher dimensional space, e.g. curves, tubular structures around curves, surfaces or 
layers between two surfaces in three dimensions, has attracted attention from many authors 
in the past three decades. Recently there has been considerable progress in rigorously, 
and non-perturbatively analyzing the changes to the spectrum of the Laplacian in curved 
geometries in two dimensions - these have been motivated by applications in mesoscopic 
semi-conductor physics, M, M. There are also studies of the contribution of torsion to the 
bound state energies, 0,0. Until now the results for higher dimensional confinement are 
restricted to the perturbative regime, f|, [|], 0, [|, f|, [1, pj . 

Below we summarize the main results of these studies, as they relate to the application 
of interest here. The next section reviews the relation between quantum mechanics and 
polymer physics in some detail. We compute the free energy in the geometries of interest 
in the third section. Most of the details of these calculations are left to an appendix. The 
final section computes the effective bending rigidities of a dilute polymer solution confined 
in a box with elastic walls, and includes our concluding remarks. 

The correct quantization procedure for quantum particles confined to manifolds embed- 
ded in higher dimensional spaces turns out to be non-trivial, [Q. Studies, || - [T2[], which are 



perturbative in the curvature of the manifold, conclude that confinement to a curved mani- 
fold results in an effective potential, which depends both on the manifold to which the particle 
is confined and the particular way this manifold is embedded in the higher dimensional space, 
i.e. both on the intrinsic and extrinsic curvatures of the manifold. This dependence on the 
extrinsic curvature means that one cannot arrive at the correct Schrodinger's equation by 
just transforming the Laplacian to the coordinates of the confining manifold; this can be 
understood in the light of the uncertainty principle - the particle will know that there is an 
"outside" , no matter how strongly it is confined. 

Confinement to a curve in three dimensions generates a potential, which is attractive to 
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the most highly curved regions, in the one dimensional Schrodinger's equation; this potential 
is proportional to the square of the local curvature of the curve, |Tl[ . If one considers a 'thin- 
tube' that is wound around a curve and has constant cross-section, (see Fig.l), to lowest order 
in the curvature one finds the same result: the energy change is proportional to the square of 
the curvature of the curve that the tube follows. The proportionality constant is independent 
of the shape of the cross-section, and also of the exact details of the confining potential (hard- 
walls etc.), 0,0, and [[U| - this can be easily understood by dimensional analysis, scaling 



arguments and using the result of ||11|| . Confinement between two equidistant surfaces in 
general dimension results in a potential U, which (in the appropriate units, i.e. neglecting 
overall positive constants) is given to second order in the principal radii of curvatures Ri, 
(z = l,2,3,.,(D-l)),by §: 

For a spherical shell embedded in D dimensions one sets Ri = a, for all i, where a is the 
radius of the shell, and, as first written down in [IIJ, Eq.(|l]) reduces to U = ( ' D ~ 1 ^2 D ~ 3 ' > ■ This 
is negative for 1 < D < 3. 

Of course, all of the above results apply to any physical system that contains the Lapla- 
cian in its equations of motion. 

II. CONNECTION TO POLYMER PHYSICS 

The particle-in-a-box problem of quantum mechanics is closely related to confined ideal 
polymers. A flexible polymer in solution is a random walk. In the case without excluded 
volume interactions this will be denoted as an 'ideal chain'; i.e. ideal chains are allowed 
to self-intersect without raising their energy. An ideal polymer chain configuration may be 
mapped onto the motion of a quantum particle, where the probability distribution function 
is the analogue of the wave function, and the degree of polymerization (N) plays the role of 
imaginary time [|15|| . In this spirit we reduce the 'ideal-chain-in-a-box' problem to a quantum 
'particle-in-a-box' problem. 



We consider a random walk in space, and work in the continuum description. Let 
G(f,f',N) be the probability for a walk of iV steps to have its ends at r and f'. The 



function G is defined to be for N < 0, and it satisfies |T5| the differential equation: 

[~M " ? V ") ^""'^ = ( 2 ) 

Here Z is the monomer length. Eq. (fj) describes diffusion once the boundary conditions are 
specified. We may also think about it as Schrodinger's equation, where N plays the role of 
imaginary time. 

The partition function is just a sum over all configurations: 

Z = (1/V) f drdr'G{r,r',N). (3) 



The factor of 1/V keeps the volume fixed when the free energy is computed. 

Here we comment on the exact problem that we are considering, and our choice of 
boundary conditions, when we compare the flat and curved boxes. Our main aim is to find 
out the effects of polymer confinement on the bending elasticities of the confining structure. 
Hence we would like to isolate the change in the free energy due to curvature only. Therefore 
we pick periodic boundary conditions in the transverse directions of the flat box. Also to 
the same end we conserve the volume of the box during bending; this, of course, will fix the 
radius of the curved box once the amount of material (the surface area of the flat box) is 
fixed. Assuming that the free energy per unit area for this system can be expanded in the 
curvature when this radius is large compared to the thickness of the box, one can compute 
the coefficients in this expansion to second order in the curvature. There are in general two 
bending elasticities to compute. These can be uniquely determined from the two equations 
one attains by considering the free energy of the sphere and the cylinder. We have in mind 
membranes which generally do not have free edges but end at the walls of a container, or 
are closed vesicles. Then, e.g. for a container with repulsive walls, one has to take into 
account when computing the energy balance the fact that the eigenfunctions for the flat box 
vanish at the transverse edges. Also, the volume for the real system may not be conserved 



- this would bring translational entropy pieces into the energy balance in addition to the 
contributions we calculate in the present work. 

To compute the partition function we use the method of eigenf unctions and levels of the 
Schrodinger's equation. This is an 'effective Hamiltonian' method where one integrates out 
certain configurations of the polymer chain and labels the coarse-grained modes by discrete 
'levels'. Thus we write Z = Y^e fc s T , where the sum is over the levels and e n 's are the 

n 

eigenvalues of the effective Hamiltonian that describes the system in terms of these levels 
only. For the remainder of this paper, in order to simplify the notation, we measure all 
energies in the units of iV/c^T, and all lengths in the units of (l/y/E), (except in the final 
results) . It becomes clear what these levels are once one employs the eigenfunction expansion 
16| for the propagator of the chain: G(r,f',N) = J2%(r)^n{r')e- NEn e(N). Here Q(N) 



is the step function assuring that G = for N < 0, and the \P n (r)'s are normalized solutions 
to V 2 \I/(r) = —E^*{f). It is easy to verify that this expansion indeed solves Eq.(|2|), (recall: 
-t^Q(N) = 5(N)). The boundary conditions are \I/(r*) = at the walls which are at \r\ — a 
and \r\ = a + d. This is a slab of thickness d in one dimension, and a spherical shell of radius 
a and thickness d in any dimension D. 

Using this eigenfunction expansion in the partition function expression, Eq.(^), we find: 
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Z = (l/V)J2 e ~ NEn I dr^* n (r) f dr'^ n (r') = (1/V) ^V^" f df^ n (r) 



(4) 



We see that in the geometry of interest all states except the spherically symmetric ones 
disappear from the sum in Eq. (|]), even though they do contribute to the propagator. (In 
D — 1, the sinusoidal solutions with an even number of peaks disappear). Hence: 

oo „ 2 

Z = (1/V) e- NE " M / drA D -^R n {r) , (5) 

n=l 

where n labels the distinct s-wave levels, R n {r) is the normalized radial wavefunction, and 
we explicitly indicate the dependence of the energies on a and d. Defining: 

Jdrr( D -VR n (r) 2 

B n (a,d) = -3 (6) 

1 ^ ^((a + d) D -a D ) K } 
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where Sd-i is the surface area of the unit sphere in D dimensions, we can write for the free 
energy: 



F = -k B TMlog ]T e- {E " Mm2/6) B n (a, d) 



, n=l 



where M is the number of chains (each containing N links) in the box. 
We now proceed with the analysis of the above expression. 

III. FREE ENERGY FOR POLYMERS IN BOXES 



(?) 



We consider a spherical annulus of large radius and small thickness in arbitrary dimension 
D. The walls are taken to be repulsive for the polymer. Physical interest lies in D = 2 
(Fig. 2), which is a cylindrical shell in three dimensions^, and D = 3 (Fig. 3), a spherical 
shell in three dimensions. 

The Laplacian in D-dimensions is V 2 = -tt + 4- — -4- Here L is the angular 

^ dr r dr r z 

momentum operator in D dimensions. After standard manipulations we find the following 
quantization condition in terms of dimensionless variables: 



Jf 



/ 1 
1 + - 

e 



(E) 



-J, 



-,i / 1 
1 + - 

e 



(E) 



e 



0. 



(8) 



where E = Ed 2 , £ 2 = v 2 + (-^^) 2 , e = (-), and v 2 is an eigenvalue of the| SO(D) symmetric 
Casimir operator, and the functions J and N are the standard Bessel functions. As explained 
earlier, spherical symmetry of the geometries considered causes all but the v — states to 
drop out of the partition function. Defining fi = 4£ 2 , (for v — 0), // = (D — 2) 2 , and labelling 
these s-wave states by n = 1,2,..., one can now proceed with a perturbative analysis by 
expanding Eq.(|8]) in powers of e = -, which results in: 



The length along the axis of the cylinder is taken to be large compared to the other dimensions, 
and therefore does not play an essential role here. 
2 SO(D) is the rotation group in D dimensions. 
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e) 2 -o 3+ 0V^)+-]- » 

This is the same perturbative result cited above. We see that the curvature lowers the 
spherically symmetric energy spectrum for 1 < D < 3, does not affect it when D — 1 
(trivial), or D = 3, and increases it everywhere else. Notice that the perturbation series 
stops for D = 1, and D = 3. Also notice that the perturbative effect is the same for all 
modes up to fourth order, this simplifies treating the full free energy of the polymer, Eq. ([?]). 
This treatment an expression for B n (a,d), the computation of which is difficult in general, 
but we were able to find an exact expression in D = 3, and a perturbative one in D = 2. 
We now summarize our results. Details are left to the appendix. 



d 2 E n (a, d) = (nrr) 2 + { - - ~ 3) 



A. Three Dimensions 

Defining Z as the partition function of the flat box, and k = (Nl 2 n 2 /Qd 2 ) we find the 
exact result for the partition function: 

Z/Z = l+g(e)(l-3f(k)) (10) 

where the functions / and g are defined in the appendix. One can expand the above 
expressions in two limits: when the typical polymer dimension, R g = yNl, is smaller or 
larger than the box width, d. The latter is the so-called ground state dominance regime 
(G.S.D.) - here k > 1; the former we will denote by "the colloid limit" (C.L.), since in this 
limit the coils look like small balls - here k < 1. In G.S.D., the corrections due to B n are 
negligible, and we find that the free energy does not have any corrections due to curvature. 
In CL. we find: 

AF = k B TM±(iY (i/f-MS- ^(Sf" + ' ' ') 

B. Two Dimensions 

In this case we find (to lowest order in e): 
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Z/Z, = - |g (l - /(*) - ^ + • ■ j (11) 

where / is the same function as for three dimensions, and h is defined in the appendix. This, 

in G.S.D. leads to: 

AF=-k B TMN±£. 

In C.L. there are no terms of order y/k, and two terms of order k which cancel each other 
exactly, and the expansion starts at third order, (this is similar to what happens for small 
hard spheres, i.e. colloids, - there one finds that there is no shift in the free energy in two 
dimensions, and a positive shift in three dimensions, [|17|]): 
AF=-k B TM^{ff{i£)W + ... 

In the strong confinement (G.S.D.) regime it does not matter to which chain the monomers 
belong - the shift in the energy is proportional to the total number of monomers, MN; in 
the C.L. regime we have an ideal gas of M objects. 

We comment on some issues before computing effective bending rigidities using the above 
results: In the non-perturbative regime we solved the eigenvalue equation, Eq. @, numeri- 
cally for the ground states in D = 2,4, and 5, for 2 < ^ < 50. The large 5 behavior of the 
solutions, (see the curves plotted in Fig. 4), are exactlyQas predicted by Eq. (|9]). Notice that 
they are all converging to the flat case value of tt 2 for very small curvature. As the curvature 
increases the behavior changes its functional form, and becomes much sharper than just a 
quadratic deviation. One can follow these curves farther towards infinite curvature, and see 
that they tend to a constant value that depends on the dimensionality. This pattern is also 
valid for n ^ 1, v = (s-wave), states. 

For surfaces in three dimensions Eq.(|l]) reduces to: U = — — j^) 2 ■ The effect of 
this potential, even after integration over the entire surface, is always attractive except for 
in two cases where it has no effect at all, the flat surface, and the sphere. Notice also that 
in D = 2 one always gets an attractive potential if the confining curves are not perfectly 



3 This is not obvious from the plot given here. 
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flat. Thus in two dimensions, -(physically we think about generalized cylinders living in 
three dimensions), the free energy is reduced by any curvature, regardless of the shape of 
the boundaries. Although we are not considering large curvatures here, it is worthwhile to 
mention that the literature cited in the introduction proves that in one and two dimensional 
structures, (tubes and strips), the ground state energy of the Laplace operator decreases 
due to curvature, and there is a bound state localized in the regions of curvature. This 
is a non-perturbative statement, and it makes very weak assumptions about the detailed 
geometry. This is important because it implies that in these systems, (at least in G.S.D.), 
the polymer would localize in regions where there is curvature. One explicit example of this 
is the numerical solution of Schrodinger's equation in oval shaped rings Jl4| . 



IV. EFFECT ON « AND k 

We conclude by calculating the contribution of this effect to the effective bending cur- 
vature elasticity of an ideal polymer confined to a fluid bilayer membrane. Let f c be the 
curvature free energy per area. Up to quadratic order in the principal radii and f c 
can be written in terms of the mean and Gaussian curvatures of the surface. In terms of 
the principal radii, the mean curvature is H = f (si + 7^)> an< ^ ^he Gaussian curvature is 
K = (4-)(4-)- A general form of f c to this order is given by fL8| : f c = 2k(H — c ) 2 + kK, 



where k and k are constants usually denoted by 'bending rigidities'. The spontaneous cur- 
vature of the surface, Co, is taken as zero for the rest of this discussion, because the film is 
assumed symmetric with respect to its mid-plane. 

To find out the effects of polymer confinement on k, and k, we first consider a vesicle 
(a spherical shell in (D = 3)), with a small curvature, i.e. a shell of radius a, with a finite 
thickness d, where d <C a. For a sphere both principal radii are equal to a, and we find 
fc = \ K (^) 2 + ^(^) 2 - Next consider a cylindrical shell (D = 2); now R\ = a, but R2 = 00. 
Thus, f c — Defining $ A f = M/V = Mj (Sd) + 0(e), S being the surface area, and using 
the results for AF for the sphere and the cylinder we conclude that: 
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5k = (A; B T)($ A/ d 3 ) Q (l - /(*) - J^)) - 
5k = (k B T){$ M d 3 ) (-± + -^(A;)) , 



2tt 2 



(12) 



(13) 



where 5k and 5k are the changes in the bending rigidities. These expressions can be expanded 
in G.S.D. and C.L., to give: 



5k = -{k B T){$> M l 



3\ J 12 W > 



Ni 



i 



9 V67r 



+ ■ ■ 



fc > 1 
fc < 1 



(14) 



5k = +{k B T){^ M l s 



J£ (f ) 2 -^^(f)-ATli ((37T - 8)*/5 



; > 1 



£) + ■■■ ; fc<1 



(15) 



We have not considered any interactions between monomers. We expect that interactions 
between chains in a dense solution would screen out the effect of the walls on the concen- 
tration profile over a very short distance and hence weaken this curvature effect. We are 
currently studying this situation. But our results may hold qualitatively for concentrations 
below the overlap concentration, To be consistent with our analysis we use Gaussian 



exponents to estimate Q* m : In G.S.D. we require NQ* m = 
C.L. the walks are three dimensional and we require NQ* m 
This leads to: 



N 



N 
R„ 



. Thus $ 



l 

Nl 2 d' 
1 



In 



m jy(3/2)j3 ■ 



5 k 



(k B T) 



12 



; k > 1 
0.026 ; k < 1 



(16) 



5k 



(k B T) 



+ < 



| ; k>i 

0.06 4 (if \k<\ 



N W - 



(17) 



Edwards theory assumes slow variations of the propagator [[II]], and is therefore valid 
over large length scales compared to the lattice spacing, I, and moderately large number of 
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steps, N. (Aside: We also need d <C a, in order to use our perturbative result). Notice that 
for k < 1 we have jj (j) 2 > 1 - which can be achieved e.g. with N ~ 100, and d/l ~ 20, then 
this term gives a factor of 4. Estimating the bare rigidities of the box to be of the order of 
k B T, we see that the effect is about 1 % for re, and 5 % for re. 

We emphasize that the major result of this article is that polymer confinement might 
reduce re, and thereby induce spontaneous curvature in the system, such as transitions from 
lamellar to bicontinuous phases or tubules etc, provided that the effect persists after the 
inclusion of self-avoidance, i.e. monomer-monomer interactions, in the analysis. We expect 
that this would be the case below the overlap threshold, but future work should address 
this issue in detail. We note here that, since there are two length scales in the problem, one 
cannot solve the self-avoiding walk problem by simple scaling arguments. 

After one understands the role of interactions, one should also consider the case of a 
polymer with a finite persistence length. We expect that in the regime where the radius 
of curvature is much larger than the persistence length the result would approach the one 
for the totally flexible polymer case we treated in the present article (modulo interactions). 
But when the two lengths are comparable, the aversity to bending of rigid polymers should 
modify the effect computed here. 

Note added in proof: As we were preparing to publish we received a preprint by E. 
Eisenriegler, A. Hanke, and S. Dietrich, in which the authors also compute bending rigidities 
of systems with chains near a repulsive surface. 
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APPENDIX A: DETAILS OF THE FREE ENERGY COMPUTATION: 



We employ Eq. (^) in three dimensions, which involves integrals that are trivial to do, 
and we find for B n the following expression: 

1 + q(e) ; n = odd 
B n (a,d) = ^{ (Al) 

-og(e) ; n = even 



n 71 



where g(e) = — \ jr^&zi - Using this in Eq. (|]) results in Eq. ([H]), where / is defined as 
follows: 

We) e~ fc " 2 

f{k) = „ f 2 (A2) 



where the symbol denotes a sum over all positive even n, (n = 2,4,6, ...); and is 
over all positive odd n, (n = 1,3,5,...). As emphasized in the text Eq. (|10|) is exact with 
the function / defined as above. This function can be expanded for large arguments (g.s.d. 
regime) as: 

/(*) = ^" 3fc (l + 0(e~ 8k )) (A3) 
For small arguments one can also expand / with a little more work. First we define the 

oo 

function A(t) = ^ e~ tn . This function has the following property: 

n=— oo 

VT . ,7T 2 



A(t) = \j- t A{ T ) (A4) 

Now we define E(k) = E (e) and 0{k) = E (o) Then f(k) = f[|. (Aside: 

-% 0{k) is the partition function for the flat strip). One can easily verify that E(Q) = |j, 
and O(0) = Also, $fj& = -\{A{Ak) - 1), and ^ = -\{A{k) - A{Ak)). Using 
Eq . (\KA}) , and the fact that A(t) = 1 + 0(e _t ), when k is large, one finds that Jq dtA(t) = 

2 

2 f dtA(4t) = 2\l~Ttk + O(e^ir), when k is small. This allows an expansion of /: 

.... 1 1 - 12tt" 2 (y/iTk- k) . _i, ... 

f(k) = - \ + 0(e k). (A5) 

o 1 — An z yixk 

This function crosses over to j e~ 3k at k ~ 0.3. 
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Notice that these expressions satisfy the consistency checks: 
e = =>■ B n (d) = -is, if n is odd, and B n = 0, if n is even. Also, k = AF = 0. 

In two dimensions the integrals are also exactly doable but the final expression for B n 
is now much more complicated - we do not bother including this here. It can, however, be 
expanded in e, to yield: 



. 1 — A e 2 + „ 2 2 e 2 + • • • ; n = odd 
B n (a,d) = -^\ 16 W (A6) 

e 2 + ■ ■ • ; n = even 

Now, when one is not in G.S.D., one in principle needs to worry about the terms that 
depend on n in the expansion of the energy spectrum, Eq.(|9]). But when the argument of 
the terms in question is larger than one, we are deep into the G.S.D. regime, whereas when 
it is smaller than one these terms can be expanded in e, giving only e 4 corrections. Thus, 
they can in practice be ignored when computing bending rigidities. Defining: 



E 



(o) 



„2 



2- 1 n 2 

we find for the partition function the result in Eq. (0). Noticing that — dh ^ = 0{k), and 

2 

h(0) = one can expand h using the expressions derived above: 

i / ■* \ vr 2 1 - 12 7r- 2 A; + 32 7r-! fci ^. _ i. . 
ft (A; = — 3 7= + O e k • A8 

Again, these expressions satisfy the above-mentioned consistency checks. These expan- 
sions for / and ft are in excellent agreement with numerical results for < k < -. The 
results for the bending rigidities follow from these expansions. 
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FIGURE CAPTIONS 



Fig.l : Twisting tube example: solid cylinder < > torus. 

Fig. 2 : D = 2 annulus. 
Fig. 3 : D = 3 annulus. 

Fig. 4 : Numerical solution for E(l/e), in two, four, and five dimensions. 
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